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Abstract. An exactly solvable position-dependent mass Schrodinger equation in two di- 
mensions, depicting a particle moving in a semi-infinite layer, is rc-cxamincd in the light of 
recent theories describing superintegrable two-dimensional systems with integrals of motion 
that are quadratic functions of the momenta. To get the energy spectrum a quadratic algebra 
approach is used together with a realization in terms of deformed parafermionic oscillator 

• operators. In this process, the importance of supplementing algebraic considerations with 
a proper treatment of boundary conditions for selecting physical wavefunctions is stressed. 

• Some new results for matrix elements are derived. This example emphasizes the interest of 
a quadratic algebra approach to position-dependent mass Schrodinger equations. 
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^ 1 1 Introduction 

Quantum mechanical systems with a position-dependent (effective) mass (PDM) have attracted 
| a lot of attention and inspired intense research activities during recent years. They are indeed 

very useful in the study of many physical problems, such as electronic properties of semiconduc- 
tors p] and quantum dots [2], nuclei [3], quantum liquids [4], 3 He clusters [5], metal clusters [6], 
^ ; etc. 

$_i ■ Looking for exact solutions of the Schrodinger equation with a PDM has become an interesting 

research topic because such solutions may provide a conceptual understanding of some physical 
phenomena, as well as a testing ground for some approximation schemes (for a list of references 
see, e.g., [7]). For such a purpose, use has been made of methods known in the constant-mass 
case and extended to a PDM context, such as point canonical transformations [El El [10], Lie 
algebraic methods |11[ 112} [T3 l I14j. as well as supersymmetric quantum mechanical (SUSYQM) 
and shape- invariance techniques [15} I16j. 

Although mostly one-dimensional equations have been considered up to now, several works 
have recently paid attention to d-dimensional problems [7J |T71 [181 021 US EH E2] • In [7J (hence- 
forth referred to as I and whose equations will be quoted by their number preceded by I), we 
have analyzed <i-dimensional PDM Schrodinger equations in the framework of first-order in- 
tertwining operators and shown that with a pair (H,Hi) of intertwined Hamiltonians we can 
associate another pair (R,Ri) of second-order partial differential operators related to the same 
intertwining operator and such that H (resp. H\) commutes with R (resp. R\). In the context 
of SUSYQM based on an sl(l/l) superalgebra, R and R\ can be interpreted as SUSY partners, 
while H and Hi are related to the Casimir operator of a larger gl(l/l) superalgebra. 
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In the same work, we have also applied our general theory to an explicit example, depicting 
a particle moving in a two-dimensional semi-infinite layer. This model may be of interest in the 
study of quantum wires with an abrupt termination in an environment that can be modelled 
by a dependence of the carrier effective mass on the position. It illustrates the influence of 
a uniformity breaking in a quantum channel on the production of bound states, as it was 
previously observed in the case of a quantum dot or a bend [23} 124] . 

From a theoretical viewpoint, our model has proved interesting too because it is solvable 
in two different ways: by separation of variables in the corresponding Schrodinger equation 
or employing SUSYQM and shape-invariance techniques. The former method relies upon the 
existence of an integral of motion L, while, as above-mentioned, the latter is based on the use 
of R. In other words, the three second-order partial differential operators H, L and R form a set 
of algebraically independent integrals of motion, which means that the system is sup erintegr able. 

Let us recall that in classical mechanics [25], an integrable system on a d-dimensional mani- 
fold is a system which has d functionally independent (globally defined) integrals of motion in 
involution (including the Hamiltonian). Any system with more that d functionally independent 
integrals of motion is called superintegrable. It is maximally superintegrable if it admits the 
maximum number 2d — 1 of integrals of motion. The latter form a complete set so that any 
other integral of motion can be expressed in terms of them. In particular, the Poisson bracket 
of any two basic integrals, being again a constant of motion, can be written as a (in general) 
nonlinear function of them. Such results can be extended to quantum mechanics [26], so that 
for quantum counterparts of maximally superintegrable systems we get (in general) nonlinear 
associative algebras of algebraically independent observables, all of them commuting with H. 

The simplest case corresponds to the class of two-dimensional superintegrable systems with 
integrals of motion that are linear and quadratic functions of the momenta. The study and 
classification of such systems, dating back to the 19th century and revived in the 1960ties [271 
[28l[29], have recently been the subject of intense research activities and substantial progress has 
been made in this area (see |^Ql SIl 1321 IHBl SQl I^Sl l^^l SSI I^Sl S3 SS| 

and references quoted therein). In particular, it has been shown that their integrals of motion 
generate a quadratic Poisson algebra (in the classical case) or a quadratic associative algebra 
(in the quantum one) with a Casimir of sixth degree in the momenta and the general form 
of these algebras has been uncovered [361 EES"! [4*6l [4"7] 148] . Algebras of this kind have many 
similarities to the quadratic Racah algebra QR(3) (a special case of the quadratic Askey-Wilson 
algebra QAW(3)) |3H 132] . They actually coincide with QR(3) whenever one of their parameters 
vanishes. The eigenvalues and eigenfunctions of the superintegrable system Hamiltonian can be 
found from the finite-dimensional irreducible representations of these algebras. The latter can 
be determined by a ladder-operator method [HU [32j [33l [34] or through a realization [35j [36] 
in terms of (generalized) deformed parafermionic operators [49] . which are a finite-dimensional 
version of deformed oscillator operators [50J. 

Since our two-dimensional PDM model belongs to this class of superintegrable systems, it is 
interesting to analyze it in the light of such topical and innovative theories. This is one of the 
purposes of the present paper, which will therefore provide us with a third method for solving the 
PDM Schrodinger equation. In such a process, we will insist on the necessity of supplementing 
algebraic calculations with a proper treatment of the wavefunction boundary conditions imposed 
by the physics of the problem - a point that is not always highlighted enough. 

Another purpose of this work is to stress the interest of a quadratic algebra approach to 
PDM Schrodinger equations. If the presence of such an algebra was already noted before in 
a one-dimensional example [51] , this is indeed - as far as the author knows - the first case where 
an algebra of this kind is used as a tool for solving a physical problem in a PDM context. 

This paper is organized as follows. In Section 2, the two-dimensional PDM model of I is 
briefly reviewed and some important comments on its mathematical structure are made in 
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conjunction with the physics of the problem. In Section 3, a quadratic algebra associated 
with such a model is then introduced and its classical limit is obtained. The finite-dimensional 
irreducible representations of the algebra are determined in Section 4. Finally, Section 5 contains 
the conclusion. 



H(k) = -d*T77^ d * ~ d vT7T^ d v + Kff^), (2-3) 



2 Exactly solvable and superintegrable PDM model 
in a two-dimensional semi-infinite layer 

In I we considered a particle moving in a two-dimensional semi-infinite layer of width ir/q, 
parallel to the x-axis and with impenetrable barriers at the boundaries. The variables x, y vary 
in the domain 

D: 0<x<oo, ~ <y<^-, 

2q 2q 

and the wavefunctions have to satisfy the conditions 

i>(0,y)=0, 4i,±lj=0. (2.1) 

The mass of the particle is m(x) = moM(x), where the dimensionless function M(x) is given by 
M(x) = sech 2 qx. (2.2) 
In units wherein h = 2niQ = 1, the Hamiltonian of the model can be written as 

M(x) x y M(x) 
where we adopt the general form (12.2) and 

V^g\x) = —q 2 cosh 2 qx + q 2 k{k — 1) csch 2 qx (2-4) 

is an effective potential. This function includes some terms depending on the ambiguity para- 
meters [52], which allow any ordering of the noncommutating momentum and PDM operators 
(see equation (12.3)). In (|2.4p . the constant k is assumed positive and we have set an irrelevant 
additive constant vq to zero. 

As shown in I, both the operators 

L = -dl 

and 

R (k) = 7? (fc)t T? (fc) 

= — cosh 2 qx sin 2 qyd 2 + 2 sinh qx cosh qx sin qy cos qy d 2 y — sinh 2 qx cos 2 qy d 2 
+ q sinh qx cosh qx(l — 4 sin 2 qy)d x + q(l + 4 sinh 2 qx) sin qy cos qyd y 
+ g ,2 (sinh 2 qx — sin 2 qy — 3 sinh 2 qx sin 2 qy) — q 2 k(l + csch 2 qx sin 2 qy) 
+ q 2 k 2 csch 2 qx sin 2 qy, 

where 

j/fc)t = — cosh qx sin qy d x + sinh qx cos qyd y — q sinh qx sin qy — qk csch qx sin qy, 
r/ fc ) = cosh qx sin qy d x — sinh cos qyd y + q sinh sin qy — qk csch sin qy, 

commute with H^ k \ although not with one another. Hence one may diagonalize either 
and L or and simultaneously. This leads to two alternative bases for the Hamiltonian 
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eigenf unctions, corresponding to the eigenvalues 

E$ = q 2 (N + 2)(N + 2k + l), N = 0,1,2,..., (2.5) 

with degeneracies 
~N 



deg(iV) 



+ 1, (2.6) 



where [N/2] stands for the integer part of N/2. 

The first basis is obtained by separating the variables x, y in the PDM Schrodinger equation 
and its members, associated with the eigenvalues (I + l) 2 q 2 of L, read 

$}(x,y) = <$}(x)xi(v), n,l = 0,1,2,..., (2.7) 
with N = 2n + I, 

^J( x ) = Af^(tanhqx) k {sechqx) l+2 pi k ~^ l+1 \l - 2tanhV), 

I— cos[(Z + l)qy] for 1 = 0,2,4,..., 
xi(y) = { y * (2.8) 
' — sin[(Z + l)qy] for I = 1, 3, 5, . . ., 

7T 

and ftf^i a normalization constant given in equation (13.18). 
The second basis, resulting from the intertwining relation 

v (k) H (k) = i? W r? (fc) ) H W = H (k+i) + 2q 2^ 
and its Hermitian conjugate, can be built by successive applications of operators of type r]( k H , 
y) =* { N %v m V {k+1) ' ■ ■■V {k+u - 1) ^t,N ^yl (2-9) 

on functions ^>^^ o (x,y), annihilated by r/ fc+z/ ) and given in Eqs. (13.28), (13.32) and (13.34). 
In (|2.9p . iVo runs over 0, 2, 4,. . . , N or N — 1, according to whether N is even or odd, while v, 
defined by v = N — Nq, determines the eigenvalue 

= q 2 u(v + 2k), u = 0,1,2,.... (2.10) 

Although an explicit expression of the normalization coefficient M^ No is easily obtained (see 

equation (13.41)), this is not the case for No (x, y) (except for some low values of N and Nq), 

nor for the expansion of ^^ No (x,y) into the first basis eigenfunctions ip^/j(x,y), which is given 
by rather awkward formulas (see equations (13.46), (13.51), (13.55) and (13.56)). 

Before proceeding to a quadratic algebra approach to the problem in Section 3, it is worth 
making a few valuable observations, which were not included in I. 

Mathematically speaking, the separable Schrodinger equation of our model admits four lin- 
early independent solutions obtained by combining the two independent solutions of the second- 
order differential equation in x with those of the second-order differential equation in y. Among 
those four functions, only the combination ip^/J(x, y), considered in (|2.7p . satisfies all the bound- 
ary conditions and is normalizable on D. It is indeed clear that the alternative solution to the 
differential equation in x is not normalizable, while that to the differential equation in y, 

_ fwx[(l + l)qy] for 1 = 0,2,4,..., 

I cos[(Z + l)qy] for I = — 1, 1, 3, 5, . . ., 
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violates the second condition in equation (|2.ip . Hence the three remaining combinations provide 
unphysical functions. 

Some mathematical considerations might also lead to another choice than L and for the 
basic integrals of motion complementing H^ k \ First of all, instead of L, one might select the 
operator p y = —id y , which obviously satisfies the condition [H^ k \p y ] = 0. This would result 
in a linear and a quadratic (in the momenta) integrals of motion, generating a much simpler 
quadratic algebra than that to be considered in Section 3. It should be realized, however, that 
the eigenfunctions e imy (m S Z) of p y , being linear combinations of the physical and unphysical 
functions (|2.8p and (12. lip , are useless from a physical viewpoint. We are therefore forced to 
consider the second-order operator L instead oip y . 

Furthermore, it is straightforward to see that another pair of first-order differential operators 

^(fe)t — _ coshgx cos qy d x — sinhqx sin qyd y — q sinhqx cos qy — qk cschqx cos qy, (2.12) 
fjW = cosh qx cos qy d x + sinh qx sin qy d y + q sinh qx cos qy — qk csch. qx cos qy , (2-13) 

intertwines with and h[ , i.e., satisfies the relation 

n {k) H {k) = jyWjjM H W = H (k+i) + 2q 2 K (2 14) 

and its Hermitian conjugate. Such operators correspond to the choice a = c = g = 0, b = d = 1 
in equation (12.29). 

As a consequence of (I2.14p . the operator 

R (k) = n (k)i-(k) 

= — cosh 2 qx cos 2 qy 9 2 — 2 sinh qx cosh qx sin qy cos qy d xy — sinh 2 qx sin 2 qy d y 
+ q sinhgx cosh(7x(l — 4 cos 2 qy)d x — q(l + 4 sinh 2 qx) sin qy cos qyd y 
+ g 2 (sinh 2 qx — cos 2 qy — 3 sinh 2 cos 2 qy) — q 2 k(l + csch 2 qx cos 2 qy) 
+ q 2 k 2 csch 2 qx cos 2 qy, 

commutes with and is therefore another integral of motion. It can of course be expressed 
in terms of H^ k \ L and R^ k \ as it can be checked that 

#( fc ) =L + R {k) + R (k) + 2q 2 k. 

However, we have now at our disposal three (dependent) integrals of motion L, anc l R( k ) 
in addition to H^ k \ so that we may ask the following question: what is the best choice for the 
basic integrals of motion from a physical viewpoint? 

This problem is easily settled by noting that the zero modes of fj^ k \ 

ui^(x,y) = (tanh qx) k (sech qx) s+1 (sin qy) s , 

violate the second condition in equation (12. ip for any real value of s and therefore lead to 
unphysical functions. This contrasts with what happens for the zero modes u>s (x, y) of 77^, 
given in (13.28), which are physical functions for s > and can therefore be used to build the 
functions ty^) No (x,y) considered in (|2.9p . as it was shown in (13.32). We conclude that the 
physics of the model imposes the choice of L and R^ as basic integrals of motion. 
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3 Quadratic associative algebra and its classical limit 

It has been shown |36} WI\ that for any two-dimensional quantum superintegrable system with 
integrals of motion A, B, which are second-order differential operators, one can construct a 
quadratic associative algebra generated by A, B, and their commutator C. This operator is 
not independent of A, B, but since it is a third-order differential operator, it cannot be written 
as a polynomial function of them. The general form of the quadratic algebra commutation 
relations is 

[A,B] = C, (3.1) 
[A,C]=aA 2 + 1 {A,B}+5A + eB + (, (3.2) 
[B, C] = aA 2 - jB 2 - a{A, B} + dA-5B + z. (3.3) 

Here {A, B} = AB + BA, 

5 = 5(H) =60 + 5^, e = e(H)=e + e 1 H, ( = ((H) = Co + (iH + ( 2 H 2 , 
d = d(H) = d + diH, z = z(H) = z + Z\H + z 2 H 2 , 

and a, 7, a, Si, e$, d, di, Z{ are some constants. Note that it is the Jacobi identity [A, [B, C]] = 
[B, [A,C]] that imposes some relations between coefficients in (j3.2|) and (|3.3p . 

Such a quadratic algebra closes at level 6 [U] or, in other words, it has a Casimir operator 
which is a sixth-order differential operator [36], 

K = C 2 + faA 3 - \a{A, A, B} - \*y{A, B, B} + (fa 2 + d + |a 7 ) A 2 
+ (§07 - 5) {A, B} + (§ 7 2 - e) B 2 + (|a<5 + |ae + |d 7 + 2z) A 
+ (-|ae + §7<5 - 2C) 5 + i 7 z - i a ( 
= k + hH + k 2 H 2 + k 3 H 3 , (3.4) 

where fc, are some constants and {A, B, C} = ABC + ACB + BAC + BCA + CAB + CBA. 

For our two-dimensional PDM model, described by the Hamiltonian defined in equations 
(MMZ3J, we shall take 

A = R, B = L, (3.5) 

where, for simplicity's sake, we dropped the superscript (k) because no confusion can arise 
outside the SUSYQM context. 

To determine their commutation relations, it is worth noting first that their building blocks, 
the first-order differential operators d y , rf and 77, generate another quadratic algebra together 
with the other set of intertwining operators ff , fj, given in (|2.12p and f|2. 13j) . Their commutation 
relations are indeed easily obtained as 

[9 y ,v] = IV, [9 y ,fj] = -qrj, [v,v] = Qd y , (3.6) 

[r,, fjt] = 2q 2 k(l + f), [fj, n ] ] = 2q 2 k(l + f ), [ V , fj*] = -qd y + 2q 2 k& (3.7) 

and their Hermitian conjugates. In (|3.7p . we have defined 

£ = —(2qk)~ 1 (rj + ry) = csch qx sin qy, f = —(2qk)~ 1 (fj + fj') = cschqxcosqy. 

Interestingly, d y , rj and fj (as well as d y , rp and fj') close an sl(2) subalgebra. 

From these results, it is now straightforward to show that the operator C in (|3.ip is given by 

C = q{d y ,rffj + fj^r}} 
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and that the coefficients in (|3.2[) and f|3.3[) are 

a = 7 = 8g 2 , 5 = 8q 2 [q 2 (2k - 1) - H], e = 16g 4 (fc - l)(fe + 1), 

^ = 8q 4 (/c - l)(2 Q 2 fc - i^), o = 0, d=16g 4 , z = 8q i {2q 2 k - H). (3.8) 

On inserting the latter in (|3.4f) . we obtain for the value of the Casimir operator 

K = -4g 4 [2g 2 (7A; - 6) - 3H](2q 2 k - H). 

It is worth noting that since a = in (|3.3p , we actually have here an example of quadratic Racah 
algebra QR(3) [E]. 

Before proceeding to a study of its finite-dimensional irreducible representations in Section 4, 
it is interesting to consider its classical limit. For such a purpose, since we have adopted units 
wherein H = 2mo = 1, we have first to make a change of variables and of parameters restoring 
a dependence on h (but keeping 2niQ = 1 for simplicity's sake) before letting h go to zero. 

An appropriate transformation is 

X = fix, Y = hy, P x = -ihdx, Py = -ihdy, Q = ?, K = hk. 

n 

On performing it on the Hamiltonian given in equations (|2.2p ~ (|2.4p . we obtain 

H = -h 2 (d x cosh 2 QXd x + dy cosh 2 QXd Y ) - h 2 Q 2 cosh 2 QX + Q 2 K(K - h) csch 2 QX, 

yielding the classical Hamiltonian 

H c = ]imH = cosh 2 QX(P X + Py) + Q 2 K 2 csch 2 QX. 
h— >o 

A similar procedure applied to the intertwining operators leads to 

r\ c = lim r] = i cosh QX sin QYPx — i sinh QX cos QYPy — QK csch QX sin QY, 
7y c = lim f\ = i cosh QAT cos QYPx + i sinh QX sin QYPy — QK csch cos QY, 

together with rj* = lim 77* and r7* = lim , while the operators quadratic in the momenta give 

H— »0 ?i^0 

rise to the functions 

L r = lim L = Py, 

R c = lim R = cosh 2 QX sin 2 QYPl - 2 sinh QA" cosh QA sin QY cos QYP x Py 

h— »o 

+ sinh 2 QA cos 2 QYPy + Q 2 if 2 csch 2 QA sin 2 QY, 

^ c = lim R = cosh 2 QA cos 2 QYP\ + 2 sinh QA cosh QA sin QY cos QYP x Py 

h— »o 

+ sinh 2 QA sin 2 QYP y + Q 2 K 2 csch 2 QA cos 2 QY, 

satisfying the relation 

H c = L c + R c + i? c . 

The quadratic associative algebra (|3.1|) - (|3.4|) is now changed into a quadratic Poisson alge- 
bra, whose defining relations can be determined either by taking the limit lim(i^,) _1 [0, 0'\ = 

{O c , 0' c }p or by direct calculation of the Poisson brackets {O c , 0' c }p: 
{A c , B c }p = C c , 
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{A C ,C C } P = a c A 2 c + 2 lc A c B c + 5 C A C + e c B c + Cc, 

{B c , C C } P = a c A 2 c - j c B 2 - 2a c A c B c + d c A c - 5 C B C + z c 



Here 



and 



C c = lim ^ = 2QP Y (r ) *fj c + fj* c Vc) 



a c = 7c = -8Q 2 , 5 C = 8Q 2 H C , e c = -16Q 4 K 2 , Cc = a c = d c = z c = 0. 



Such a Poisson algebra has a vanishing Casimir: 
K c = lim K = 0. 

h-,0 

4 Finite-dimensional irreducible representations 
of the quadratic associative algebra 

The quadratic algebra (I3.ip - (|3.4p can be realized in terms of (generalized) deformed oscillator 
operators J\f, b^ , b, satisfying the relations [50] 



[AT, = b\ [N, b] = -b, tfb = $(M), btf = $(AA + 1), 

where the structure function §(x) is a 'well-behaved' real function such that 

$(0) = 0, > for x > 0. (4.1) 

This deformed oscillator algebra has a Fock-type representation, whose basis states \m), m = 0, 
1,2, . . .0 fulfil the relations 

AA|m) = m\m), 

tf\m) = y/^im + 1) \m + 1), m = 0,1,2,..., 

6|0) = 0, (42) 



6|m) = y $(m) |m — 1), m = l,2,.... 

We shall be more specifically interested here in a subclass of deformed oscillator operators, 
which have a (p + l)-dimensional Fock space, spanned by \p, m) = \m), m = 0, 1, . . . ,p, due to 
the following property 

$(p + l) = (4.3) 
of the structure function, implying that 
(&t)P+i = feP+i = o. 

These are so-called (generalized) deformed parafermionic oscillator operators of order p 
The general form of their structure function is given by 

= x(p + 1 — x)(a + aix + a 2 x 2 H h ap.ix^" 1 ), 



x We adopt here the unusual notation m) in order to avoid confusion between the number of deformed bosons 
and the quantum number n introduced in (|2.7| | . 
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where do, Oi, • • • , may be any real constants such that the second condition in (|4.ip is 
satisfied for x = 1, 2, ... ,p. 

A realization of the quadratic algebra (|3.ip - (|3.4p in terms of deformed oscillator operators 
M,b\b reads [SB] 

(4.4) 
(4.5) 



A = A{M), 

B = cr(Af) + tfp{N) + p{M)b, 



where A(Af), cr(Af) and p{M) are some functions of Af, which, in the 7 ^ case, are given by 



A(Af)- 
a{M)- 
p 2 {Af) 



1 



(Af + u) 



a 
4 



(Af + u) 



+ 



ae — 7<5 ae 2 — 2^5e + 4:j 2 ( 
~~ 2^ 



4 7 4 



(jV + u) 



1 • 



1 



3 • 2 12 7 8 (AA + u)(jV + u + l)[2(N + u) + l] 2 ' 
with the structure function 



(4.6) 
(4.7) 
(4.8) 



$(x) = -3072 7 6 ^:[2(AA + u) - l] 2 

- 48 7 6 (a 2 e - a-f5 + a 7 e - (i7 2 )[2(AA + u) - 3][2(AA + u) - l] A [2{M + u) + 1] 
+ 7 8 (3a 2 + 4a 7 )[2(AA + u) - 3] 2 [2(AA + u) - \} 4 {2{M + u) + l] 2 

+ 768(ae 2 - 2j5e + 4 7 2 C) 2 

+ 32 7 4 (3a 2 e 2 - 6a-f5e + 2a7e 2 + 2 7 2 <5 2 - id^e + 87 3 z + 4a7 2 C) 
x [2(Af + u) - l] 2 [12(AA + n) 2 - 12(Af + u) - 1] 

- 2567 2 (3a 2 e 3 - 9a7<5e 2 + a^e" + 67 2 5 2 e - M^fe 2 + 2 7 4 5 2 + 2cfy 4 e + 12 7 3 e; 

- 4 7 5 z + 12a 7 2 eC - 12 7 3 5C + 4a 7 4 C)[2(AA + u) - l] 2 . 



(4.9) 



These functions depend upon two (so far undetermined) constants, u and the eigenvalue of the 
Casimir operator K (which we denote by the same symbol). 

Such a realization is convenient to determine the representations of the quadratic algebra 
in a basis wherein the generator A is diagonal together with K (or, equivalently, H) because 
the former is already diagonal with eigenvalues A{m). The (p + l)-dimensional representations, 
associated with (p + l)-fold degenerate energy levels, correspond to the restriction to deformed 
parafermionic operators of order p [36]. The first condition in (|4.ip can then be used with 
equation (|4.3|) to compute u and K (or E) in terms of p and of the Hamiltonian parameters. A 
choice is then made between the various solutions that emerge from the calculations by imposing 
the second restriction in (|4.1|) for x = 1,2, ... ,p. 

In the present case, for the set of parameters (|3.8p . the complicated structure function (|4.9p 
drastically simplifies to yield the factorized expression 



$(x) = 3 • 2 3 V°(2:e + 2u + k - l)(2x + 2u + k - 2) (2a + 2u - k)(2x + 2u - k - 1) 

x (2x + 2u - \ + A) (2x + 2u - § + A) (2x + 2u - \ - A) (2x + 2u - § - A) , 



where 
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Furthermore, the eigenvalues of the operator A become 

A(m) = q 2 (2m + 2u- k)(2m + 2u + k). 

Since A = R is a positive-definite operator, only values of u such that A(m) > for m = 0, 
1, . . . ,p should be retained. 

On taking this into account, the first condition in (|4.ip can be satisfied by choosing either 
u = k/2 or u = (k + l)/2, yielding 

A(m) = Aq 2 m{m + k) (4-10) 

or 

A(m) = 4g 2 (m + |) (m + /c + \) , (4.11) 

respectively. For u = k/2, equation (14. 3h . together with the second condition in (14. ip . can be 
fulfilled in two different ways corresponding to A = 2p + k + 1 db \ or 

£ = g 2 (2p + §±±) (2p + 2fe + i±i). (4.12) 

The resulting structure function reads 

= 3 • 2 38 g 20 a;(p + 1 - x) (s - \) (p + 1 ± | - s) (x + k - \) (x + k - 1) 

x ( x + p + k + l±i) (x+p + k-\±\). (4.13) 

Similarly, for u = (k + l)/2, we obtain 

^ = g 2 (2p+|±|) (2p + 2A; + |±i) (4.14) 

and 

$(x) = 3 • 2 38 g 20 x(p + 1 - x) (x + \) (p + 1 ± \ - x) (x + k) (x + k - \) 

x (x + p + k + | ± \) (x + p + k + | ± i) . (4.15) 

Our quadratic algebra approach has therefore provided us with a purely algebraic derivation 
of the eigenvalues of H and R in a basis wherein they are simultaneously diagonal. It now 
remains to see to which eigenvalues we can associate physical wavef unctions, i.e., normalizable 
functions satisfying equation (|2.ip . This will imply a correspondence between \p,m) and the 
functions ^jv,jv— i/(a?5 defined in (|2.9p . 

On comparing A{m) to the known (physical) eigenvalues r v of R, given in (|2.10p . we note 
that the first choice (I4.10|) for A(m) corresponds to even u = 2m (hence to even N), while the 
second choice (|4.1ip is associated with odd v = 2m + 1 (hence with odd N). Appropriate values 
of p leading to the level degeneracies (12. 6h are therefore p = N/2 and p = (N — l)/2, respectively. 
With this identification, both equations (|4.12p and (|4.14p yield the same result 

E = q 2 (N +l±\) (N + 2k + \±\) . (4.16) 

Comparison with (|2.5p shows that only the upper sign choice in (|4.16p leads to physical wave- 
functions ^N,N-u(x,y). 

Restricting ourselves to such a choice, we can now rewrite all the results obtained in this 
section in terms of iV and v instead of p and m. In particular, the two expressions (I4.13P 
and (|4.15p for the structure function can be recast in a single form $(m) — > & u , where 

§ v = 3 • 2 3 W^ - + 2k - l)(y + 2k - 2)(N + v + 2k) {N + v + 2k + 1) 

x (N- u + 2){N- i/ + 3). (4.17) 
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More importantly, our quadratic algebra analysis provides us with an entirely new result, 
namely the matrix elements of the integral of motion L in the basis wherein H and R are 
simultaneously diagonal. On using indeed the correspondence \p, m) — ► N— v, as weu as the 
results in equations (jOJ) . (j4"3|) , (j4"77]) . P~8j) and (f^7T7|) . we obtain 



(4.18) 



where we have reset a(m) — ► a u , p(m) — > p u and defined t v = s u p v ^2\f^u- The explicit form of 
the coefficients on the right-hand side of (|4.18p is given by 



-{-{v + k-l) 2 {v + k + lf 



2{y + k- \){v + k + 1) 
+ [N 2 + {2k + 3)N + 2k 2 + 2k + \}{v + k- l)(v + k + 1) 
- k(k - 1){N + k + 1)(N + k + 2)}, 

-1/(1/ - \){y + 2k - l)(v + 2k - 2) 



16(z/ + A; - 2)(z/ + fe - l) 2 (z/ + it) 
x (JV - v + 2)(iV - v + 3)(iV + i/ + 2k){N + z/ + 2k + 1). 



(4.19) 



(4.20) 



Note that t v is determined up to some phase factor s u depending on the convention adopted for 
the relative phases of ^>n,n~u and ^ n,n-u+2- 

For N = 4, for instance, v runs over 0, 2, 4, so that equations (|4.18p - (|4.20p become 



4,0 



. 2(fc + 1)(2A; + 3)(2fe + 9) _ 
(13fc + 21)^,0 + 3s 4 \/ 1 fc + 2 ^ 4 > 2 



L$ 4 2 = q 2 



k + 3 



3s 4 /2(/c + l)(2k + 3)(2/c + 9) 17/c 2 + 76/c + 39 T 

-w 4i o H — r~rT7i n\ ^4,2 



k + 2 



(k + l){k + 3) 



k + 1 



s 2 10(k + 3)(2k + l)(2k + 7) 
+ - A; r— ; V 4 a 



k + 2 



k + 1 



10(k + 3)(2A; + l)(2k + 7) D , _ 

\l — \ + 2 4,2 ^ ^ 4,4 



As a check, these results can be compared with those derived from the action of L on the 
expansions of Vl/^o, ^4,2 and ^4^ in terms of the first basis eigenfunctions ^0,4) ^1,2 and "02,0 
(see, e.g., equations (13.61) and (13.49) for ^4^ and ^4,4, respectively). This leads to the phase 
factors S2 = S4 = — 1. 

To conclude, it is worth mentioning that had we made the opposite choice in equation (|3.5p . 
i.e., A = L and B = R, we would not have been able to use the deformed parafermionic 
realization (|4.4p . (|4.5p to determine the energy spectrum. The counterpart of the parafermionic 
vacuum state would indeed have been a function annihilated by L and therefore constructed 
from the unphysical function x-i(y) of equation (|2.1ip . 



5 Conclusion 

In this paper, we have revisited the exactly solvable PDM model in a two-dimensional semi- 
infinite layer introduced in I. Here we have taken advantage of its superintegrability with two 
integrals of motion L and R that are quadratic in the momenta to propose a third method of 
solution in the line of some recent analyses of such problems. 
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We have first determined the explicit form of the quadratic associative algebra generated 
by L, R and their commutator. We have shown that it is a quadratic Racah algebra QR(3) and 
that its Casimir operator K is a second-degree polynomial in H. We have also obtained the 
quadratic Poisson algebra arising in the classical limit. 

We have then studied the finite-dimensional irreducible representations of our algebra in 
a basis wherein K (or H) and R are diagonal. For such a purpose, we have used a simple proce- 
dure, proposed in [36], consisting in mapping this basis onto deformed parafermionic oscillator 
states of order p. Among the results so obtained for the energy spectrum, we have selected 
those with which physical wavefunctions can be associated. This has illustrated once again the 
well-known fact that in quantum mechanics the physics is determined not only by algebraic 
properties of operators, but also by the boundary conditions imposed on wavefunctions. Our 
analysis has provided us with an interesting new result, not obtainable in general form in the 
SUSYQM approach of I, namely the matrix elements of L in the basis wherein H and R are 
simultaneously diagonal. 

As final points, it is worth observing that the approaches followed here are not the only 
ones available. First, one could have used a gauge transformation to relate equation (|2.3p 
to a well-known superintegrable system in a Darboux space ( |38[ HH] and references quoted 
therein). Second, the irreducible representations of QR(3) could have been constructed by the 
ladder-operator method employed in [3H [32j [331 GS] • This would have allowed us to express the 
transformation matrix elements between the bases ip^\ and ' N (denoted by Z^. n l in I) in 
terms of Racah- Wilson polynomials. 
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